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We present a search for standard model Higgs boson production in association with a W boson 
in proton-antiproton collisions [pp W^H — > £ubb) at a center of mass energy of 1.96 TeV. The 
search employs data collected with the CDF II detector that correspond to an integrated luminosity 
of approximately 1.9 fb^^. We select events consistent with a signature of a single charged lep- 
ton (e*//i*), missing transverse energy, and two jets. Jets corresponding to bottom quarks are 
identified with a secondary vertex tagging method, a jet probability tagging method, and a neural 
network filter. We use kinematic information in an artificial neural network to improve discrimina- 
tion between signal and background compared to previous analyses. The observed number of events 
and the neural network output distributions are consistent with the standard model background 
expectations, and we set 95% confidence level upper limits on the production cross section times 
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branching fraction ranging from 1.2 to 1.1 pb or 7.5 to 102 times the standard model expectation 
for Higgs boson masses from 110 to 150GeV/c^, respectively. 



PACS numbers: 13.85.Rm, 14.80.Bn 
I. INTRODUCTION 

Standard electroweak theory predicts a single funda- 
mental scalar particle, the Higgs boson, which arises as a 
result of spontaneous electroweak symmetry breaking [l[ . 
However, the Higgs boson is the only fundamental stan- 
dard model particle which has not been directly observed 
by experiments. The current experimental constraint 
on the Higgs boson mass, rriff > 114.4 GeV/c^ at 95% 
confidence level (C.L.), comes from direct Higgs boson 
searches at LEP experiments 0. Global fits to elec- 
troweak measurements exclude masses above 154 GcV/c^ 
at 95% confidence level [1]. 

At the Tevatron pp collider at Fcrmilab, the ncxt-to- 
leading-order (NLO) Higgs boson production cross sec- 
tion is about 10 times larger for gluon fusion than for 
WH associated production, and the cross section for 
WH is about twice that of ZH 0]. The Higgs boson 
decay branching fraction is dominated hy H bb for 
niH < 135GeV/c2 and by H W+W' for tuh > 
135GeV/c^ Background bb production processes 

have QCD cross sections at least 4 orders of magnitude 
greater than that of Higgs boson production @ , and this 
renders searches in the gg H bb channel nonvi- 
able. Requiring the Icptonic decay of the associated W 
boson reduces the huge bb background rate. As a result, 
WH Ivbb is one of the most favorable channels for a 
low mass Higgs boson search^ , and it forms an important 
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component of the combined search for the Higgs boson 
at the Tevatron. 

Searches for WH ivbb at ^/s = 1.96 TeV have been 
most recently reported by CDF 0, |1| and DO ^ us- 
ing data corresponding to an integrated luminosity of 
955 pb~^ and 440 pb^^, respectively. In this paper, we 
present an update on the search for WH £vbb produc- 
tion at CDF using about 1.9 fb~^ of data and improved 
analysis techniques. The b hadrons which arise from bb 
quark pairs in this final state can be identified through 
special algorithms based on the tracking information for 
particles within hadronic jets, a procedure known as b- 
tagging. We have optimized the ^-tagging algorithms to 
increase the acceptance for double ^-tagged events. In 
addition we increase the signal acceptance by including 
the electrons going into the forward region of the detec- 
tor and introduce a multivariate discriminant technique 
using a neural network (NN) to reduce large background 
contamination after the event selection. 

The paper is organized as follows. Section HIl describes 
the CDF II detector. The event selection criteria are ex- 
plained in Sec. IIIII In Sec. IIVI the 6-tagging algorithms 
SECVTX neural net 6-tagging filter, and jet proba- 
bility [ll| are discussed in detail. Contributions from the 
standard model (SM) background in the WH tvbb 
channel are calculated in Sec. |V] for various sources. In 
Sec. I VII signal acceptance and systematic uncertainties 
are estimated. The neural network discriminant tech- 
nique is described in Sec. IVIII The results and their sta- 
tistical interpretation are presented in Sec. IVIIII Finally, 
our conclusions are presented in Sec. IIXI 



II. CDF II DETECTOR 

The CDF II detector [l2| geometry is described us- 
ing a cylindrical coordinate system. The z-axis follows 
the proton direction, and the polar angle is usually ex- 
pressed through the pseudorapidity -q — — ln(tan(6'/2)). 
The detector is approximately symmetric in 7] and in the 
azimuthal angle 0. The transverse energy is defined as 
Et = E sin 9, and the transverse momentum pT = p sin 9. 

Charged particles are tracked by a system of silicon 
microstrip detectors and a large open cell drift cham- 
ber in the region \'q\ < 2.0 and |?7| < 1.0, respectively. 
The tracking detectors are immersed in a 1.4 T solcnoidal 
magnetic field aligned coaxially with the incoming beams, 
allowing measurement of charged particle momentum 
transverse to the beamline (pr)- 



{Vfi) or their antiparticles. 
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The transverse momentum resolution is measured to 
be 8pt/pt ~ 0.07% • pT(GeV) for the combined tracking 
system. The resolution on the track impact parameter 
(do), the distance from the bcamline axis to the track 
at the track's closest approach in the transverse plane, 
is cr^do) Ri 40 /im, of which about 30 /im is due to the 
transverse size of the Tevatron beam itself. 

Outside of the tracking systems and the solenoid, 
segmented calorimeters with projective tower geome- 
try are used to reconstruct electromagnetic showers and 
hadronic jets [H, [13, over the pseudorapidity range 
|?7| < 3.6. A transverse energy is measured in each 
calorimeter tower where the polar angle (9) is calculated 
using the measured z position of the event vertex and the 
tower location. 

Contiguous groups of calorimeter towers with signals 
arc identified and summed together into an energy clus- 
ter. Electron candidates are identified in the central elec- 
tromagnetic calorimeter (CEM) or in the forward, known 
as the plug, electromagnetic calorimeter (PEM) as iso- 
lated, mostly electromagnetic clusters that match a track 
in the pseudorapidity range \r]\ < 1.1 and 1.1 < |?7| < 2.0, 
respectively. The electron transverse energy is recon- 
structed from the electromagnetic cluster with a preci- 
sion ct{Et)/Et « 13.5%/V£^T(GeV) ® 2% for central 
electrons [H and a{ET)/ET = 16.0%/ ^/E^iGeV) (3 2% 
for plug electrons. Jets arc identified as a group of 
electromagnetic (EM) and hadronic (HAD) calorime- 
ter clusters which fall within a cone of radius AR w 
y/ (A0)2 -I- [Arjy < 0.4 units around a high-i?T seed 
cluster [l^. Jet energies are corrected for calorimeter 
nonlinearity, losses in the gaps between towers and mul- 
tiple primary interactions. The jet energy resolution is 
approximately a{ET) « [O.IEt + 1.0 GeV] 

Muon candidates are detected in three separate sub- 
detectors. After at least five interaction lengths in the 
calorimeter, the large angle muons first encounter four 
layers of planar drift chambers (CMU), capable of de- 
tecting muons with pt > 1.4GeV/c [l^. Four additional 
layers of planar drift chambers (CMP) behind another 
60 cm of steel detect muons with pt > 2.8 GeV/c [l9l |. 
These two systems cover the same central pseudorapidity 
region with {rjl < 0.6. Muons that exit the calorimeters 
at 0.6 < \r]\ < 1.0 are detected by the CMX system of 
four drift layers. Muon candidates are then identified 
as isolated tracks which extrapolate to line segments or 
"stubs" in the muon subdctectors. A track that is linked 
to both CMU and CMP stubs is called a CMUP muon. 

The missing transverse energy (I/t) is a reconstructed 
quantity that is defined as the absolute value of the op- 
posite of the vector sum of all calorimeter tower energy 
depositions projected on the transverse plane. It is often 
used as a measure of the sum of the transverse momenta 
of the particles that escape detection, most notably neu- 
trinos. To be more readily interpretable as such, the raw 
I^T vector is adjusted for corrected jet energies and the 
muon momentum is also added for any minimum ionizing 
high-pT muon found in the event. 



Muon and electron candidates used in this analysis are 
identified during data taking with the CDF trigger sys- 
tem, a three-level filter, with tracking information avail- 
able at the first level [20|. Events used in this analysis 
have all passed the high-energy electron or muon trigger 
selection. The first stage of the central electron trig- 
ger (CEM) requires a track with pT > 8 GeV/c point- 
ing to a tower with Et > 8 GeV and i?HAD/£^EM < 
0.125. As appropriate for selecting VF-decay electrons, 
the plug electron trigger (MET-I-PEM) requires a tower 
with Et > S GeV, E-had/Eem < 0.125 and the miss- 
ing transverse energy I/t> 15 GeV. The first stage of 
the muon trigger requires a track with > 4 GeV/c 
(CMUP) or 8 GeV/c (CMX) pointing to a muon stub. 
A complete lepton reconstruction is performed online in 
the final trigger stage, where we require Et > 18 GeV for 
central electrons (CEM), Et > 18 GeV and ^t> 20 GcV 
for plug electron (MET+PEM) and pT > 18GcV/c for 
muons (CMUP,CMX). 



III. EVENT SELECTION 

The results presented here use data collected between 
February 2002 and May 2007. The data collected using 
the CEM, CMUP and MET-f PEM triggers correspond to 
1.92 ± 0.12 fb~^ of integrated luminosity, while the data 
from the CMX trigger corresponds to 1.88 ± 0.11 fb~^ 

The observable final state from the WH — > ivbb signal 
consists of two &-jets from Higgs decay, while the decay of 
W yields the high-p^ lepton and large missing transverse 
energy from the neutrino. Therefore events are consid- 
ered as WH candidates only if they have exactly one 
lepton candidate [2l| , with Et > 20 GeV for electrons or 
Pt > 20GeV/c for muons. Because the lepton from a 
Icptonic W decay is well-isolated from the rest of event, 
the cone of AR = 0.4 surrounding the lepton is required 
to contain less than 10% of the lepton energy. A primary 
event vertex position is calculated by fitting a subset of 
particle tracks that are consistent with having come from 
the beamline. The closest distance zq along the beam line 
between the primary event vertex and the lepton track 
must be less than 5 cm to ensure the lepton and the 
jets come from the same hard interaction. Some leptonic 
Z decays would mimic the single-lepton signature if one 
of the leptons is missed. Therefore Z boson candidate 
events are rejected if a track, an EM cluster or a jet to- 
gether with the primary lepton forms an invariant mass 
between 76 and 106 GeV/c^. The selected events are re- 
quired to have I/t greater than 20 GeV, to be consistent 
with the presence of a neutrino from W decay. 

In the plug region, we also require a high-pT isolated 
lepton candidate with Et > 20 GcV, with the same se- 
lection criteria as for the central region. In addition, be- 
cause the background contamination from multijet pro- 
cesses is higher in the forward region, we impose stricter 
criteria on the missing transverse energy, which improves 
the rejection of hadronic background events by a fac- 
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tor of four while only reducing the signal efficiency by 
20%. We require that METsig > 2, ^t> 25 GeV 
and I^T> 45 GcV when the I^t is pointing close to a 
jet, and large transverse mass of the reconstructed W , 
Mt{W) > 20 GeV/c2. Here, MET^ig is defined as the 
ratio of I^t to the square root of a weighted sum of jet 
Et with factors correlated with mismeasurement, such as 
angles between the I^t and the jet and size of jet energy 
corrections [1^, and Mt{W) is defined as follows: 

Mt{W) = ^2{p'^Pl^T-I^?'-m, (1) 

The WH signal includes two jets originating from H 
bb decays; these jets are expected to have large transverse 
energy. The jets are required to be in the pseudorapidity 
range covered by the silicon detector so that secondary 
vertices from b decays can be reconstructed. Specifically, 
we require the jets to satisfy Et > 20 GeV and |7/| < 2.0 
and use loose jets {Et > 12 GeV and \r]\ < 2.4) as input 
for the neural network discriminant discussed in Sec. lVIll 
The search for WH Ivbb is performed in the sample 
of events with W+ exactly 2 jets; however, samples of 
events with H^+l,3,>4 jets are used to cross-check the 
background modeling. 

To increase the signal purity of the W^+2 jets events, at 
least one jet must be 5-tagged by the SECVTX algorithm 
as explained in the subsequent sections. Three exclusive 
6-tagged event categories arc considered. The first cate- 
gory (ST-fST) is for events where there are two SECVTX 
6-tagged jets. The second category (ST-I-JP) consists 
of events where only one of the jets is 6-tagged by the 
SECVTX and the second jet is 6-tagged by jet probabihty. 
The third category (ST with NN filter) contains events 
where only one of the jets is 6-tagged by the SECVTX and 
also passes the neural network 6-tagging filter. 

These three categories are selected exclusively and 
the tightest ST-I-ST category is preferentially selected. 
Subsequently the ST-f JP category is considered if the 
ST-hST selection failed and finally ST with NN filter 
category is considered if both double 6-tagged selections 
failed. 



IV. h JET IDENTIFICATION ALGORITHMS 

Multijet final states have dominant contributions from 
QCD light-flavor jet production, but the standard model 
Higgs boson decays predominantly to bottom quark 
pairs. Correctly identifying the 6 quark jets helps to re- 
move most of the QCD background. 

The 6-quark has a relatively long lifetime, and 6 
hadrons formed during the hadronization of the initial 
b quark can travel a significant distance before decaying 
into a collection of lighter hadrons. Jets containing a 6- 
quark decay can be reconstructed by identifying tracks 
significantly displaced from the pp interaction point (pri- 
mary vertex). 

In this analysis, we employ three 6-identification algo- 
rithms to optimize the selection of 6-quark jets. 



A. Secondary Vertex 6- Tagging 

The SECVTX 6-tagging algorithm is applied to each 
jet in the event, using only the tracks which are within 
■q-cj) distance of Ai? = 0.4 of the jet direction. Dis- 
placed tracks in jets are used for the SECVTX recon- 
struction and are distinguished by a large impact pa- 
rameter significance (|do/(T(io|), where c?o and ad„ are the 
impact parameter and the total uncertainty from track- 
ing and beam position measurements, respectively. Sec- 
ondary vertices are reconstructed with a two-pass ap- 
proach which tests for high-quality vertices in the first 
pass and allows lower-quality vertices in the second pass. 
In pass 1, at least three tracks arc required to pass loose 
selection criteria [pt > 0.5GeV/c, \do/ado\ > 2.0), and 
a secondary vertex is fit from the selected tracks. One of 
the tracks used in the reconstruction is required to have 
Pt > l.OGeV/c. If pass 1 fails, then a vertex is sought 
in pass 2 from at least two tracks satisfying tight selec- 
tion criteria {pt > l.OGeV/c, jdo/o'dol > 3.5 and one of 
the pass 2 tracks must have pT > 1.5GeV/c). If either 
pass is successful, the transverse distance {Lxy) from the 
primary vertex of the event is calculated along with the 
associated uncertainty This uncertainty ctl^ includes 
the uncertainty on the primary vertex position. Finally 
jets are tagged positively or negatively depending on the 
Lr^y significance (L^y/aL^y)- 

Lxy/(JL^y > 7.5 (positive tag) (2) 
Lxy/(JL^y < -7.5 (negative tag) (3) 

This value has been tuned for optimum efficiency and pu- 
rity in simulated 6-jet samples from decays of top quarks. 
The energy spectrum for those jets is similar to the spec- 
trum for 6 jets from decays of Higgs bosons. 

The sign of L^y indicates the position of the secondary 
vertex with respect to the primary vertex along the di- 
rection of the jet. If the angle between the jet axis and 
the vector pointing from the primary vertex to the sec- 
ondary vertex is less than 7r/2, L^y is positively defined; 
otherwise, it is negative. If L^y is positive, the secondary 
vertex points towards the direction of the jet, as in true 
6 hadron decays. For negative L^y the secondary vertex 
points away from the jet; this may happen as a result 
of mismeasured tracks. In order to reject secondary ver- 
tices due to material interaction, the algorithm vetoes 
two-track vertices found between 1.2 and 1.5 cm from 
the center of the silicon detector (the inner radius of the 
bcampipe and the outer radius of the innermost silicon 
layer being within this range). All vertices more than 2.5 
cm from the center are rejected because 6-jets in WH 
events travel 3.5 mm on average. 

The negative tags are useful for evaluating the rate 
of false positive tags, which are identified as "mistags" 
in the background estimates. Mismeasurements are ex- 
pected to occur randomly; therefore the L^y distribution 
of fake tags is expected to be symmetric around zero. 
Simulated events are used to correct a small asymmetry 
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due to true long-lived particles in light-flavor jets and to 
interactions in the detector material. 

The efficiency for identifying a secondary vertex is dif- 
ferent in simulation and data. We measure an efficiency 
correction factor, which is defined as the ratio of the ob- 
served to the simulated efficiencies, to be 0.95 ± 0.04 in 
a sample of high- i?T jets enriched in b jets by requiring 
a soft lepton (p-r > 8GeV/c) from semileptonic heavy 
quark decays |10l |. 



B. Neural Network 6- Tagging Filter 

The sample tagged by the SECVTX algorithm still has 
significant contamination from falsely-tagged light-flavor 
or gluon jets and the misidentiflcation of c quarks as b- 
jets (23j . This search uses a multivariate neural network 
^N) technique to improve the SECVTX tagging purity 0, 

The neural network used in this article employs the 
JETNET [i^l package. The tagger is designed with two 
neural networks in series. The b — I network is trained to 
separate 6-jets from light-quark jets (/-jets), and the b — c 
network is trained to separate &-jets from c-jets. Jets that 
pass a cut on both of the neural network outputs are ac- 
cepted by the tagger. These neural networks are trained 
and applied only to jets that are already tagged by the 
SECVTX algorithm. The current neural network 6-tagging 
filter is tuned to increase the purity of the SECVTX b- 
tagged jets. 

The neural networks take as input 16 variables cho- 
sen to distinguish 6-quark jets from c- and /-quark jets 
by means of their higher track multiplicity, larger invari- 
ant mass, longer lifetime, and a harder fragmentation 
function than c- and /-quark jets. The track parameters 
and Lxy significance are good discriminators for 6-jcts. 
The sum of transverse momentum and mass A/vtx of 
the tracks associated with the displaced vertex are use- 
ful variables for identifying /-jets; however c-jets have pt 
spectra similar to 6-jets. Pseudo-cr {L^y x M^t^/p^^), 
the vertex fit x^, and the track-based probability for a jet 
to originate away from the primary vertex are the best 
discriminators for 6-jcts. 

The neural network 6-tagging filter is validated by com- 
paring the performance on a 6-enriched sample of SECVTX 
tagged heavy-fiavor jets from events with an electron 
candidate with Et > 8 GeV and from the correspond- 
ing Monte Carlo sample. A good agreement is found in 
neural network ^-tagging filter performance between data 
and Monte Carlo [J |^] . 

The output of the neural network is a value ranging 
from to 1. We cut on this value so that we reject 65% of 
light-flavor jets and about 50% of the c jets while keeping 
90% of &-jets which were tagged by SECVTX. The data-to- 
Monte-Carlo correction factor, measured from the elec- 
tron sample, is 0.97±0.02. Note that this is an additional 
correction factor with respect to the SECVTX efficiency 
scale factor because all of the jets under consideration 



have already been tagged by SECVTX. 



C. Jet Probability 6-tagging 

The jet probability 6-tagging algorithm differs from 
SECVTX in that it employs the signed impact parame- 
ters (and their uncertainties) of all tracks within a jet 
and calculates the probability that the jet was produced 
at a position consistent with the primary vertex. The 
signed impact parameter of a track is defined as positive 
if the angle between the jet direction and the line joining 
the primary vertex to the point of closest approach of 
the track is less than 90 degrees, and it is defined as neg- 
ative otherwise. A feature of this algorithm is that the 
6-tagging is performed using a continuous variable in- 
stead of a discrete object like a reconstructed secondary 
vertex. 

For a light-quark jet, most particles should originate 
from the primary vertex. Because of the finite tracking 
resolution, charged particle tracks can be reconstructed 
with a nonzero impact parameter and have an equal prob- 
ability to be either positively or negatively signed. Since 
a long-lived particle will travel some distance along the 
jet direction before decaying, its decay products will pref- 
erentially have positively signed impact parameters. 

To calculate the track probability, the tracking resolu- 
tion can be extracted by fitting the negative side of the 
signed impact parameter distribution from the inclusive 
jet data which is dominated by prompt jets. Tracks are 
sorted into different categories (77, px of tracks, and qual- 
ity of silicon detector hits) to parametrize their proper- 
ties. To minimize the contribution from badly measured 
tracks with large reconstructed impact parameters, the 
distribution of a related quantity, the signed impact pa- 
rameter significance (ratio of the signed impact param- 
eter to its uncertainty) is parametrized for each track 
category. Each track is required to satisfy the quality 
criteria of px > 0.5 GeV /c and a minimum number of 
hits in the tracking detector. 

To calculate the jet probability from tracks, at least 
two tracks with positive impact parameter are required. 
By definition, the jet probability distribution should be 
flat between and 1 for jets having only prompt tracks 
and peaked near zero for heavy flavor jets having large 
positively signed impact parameter. Tracks with a neg- 
ative impact parameter are used to define a negative jet 
probability, which is used to check the algorithm and to 
estimate the misidentiflcation rate. A jet is considered as 
tagged if the jet probability is less than 5%. We choose 
this value such that we reject 95% of light-flavor jets 
while keeping 60% of 6-jets. The difference in efficiency 
between the simulated and observed data is taken into 
account as a scale factor. We measure the scale factor to 
be 0.85±0.07 in a sample of high-i?r jets enriched in b 
jets by requiring a soft lepton from semileptonic heavy- 
flavor decay. A more detailed description of the scale 
factor estimation is given in Ref . [ll| . 
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V. BACKGROUND 

The final state signature of WH ivhb production 
can be mimicked by other processes. The dominant 
background processes are VF+jets production, produc- 
tion, and non-PF QCD multijet production. Several elec- 
troweak production processes also contribute but with 
smaller rates. In the following subsections the contribu- 
tion from each background source is discussed in detail. 
These background estimations are based on the same 
strategies used in the previous analysis 0, . The quan- 
titative background estimates can be found in Sec. IVIIII 



A. Non-VK QCD Multijet 

Events from QCD multijet production sometimes 
mimic the W^-boson signature due to fake leptons and 
fake I^T- A non-VF lepton is reconstructed when a jet 
passes the lepton selection criteria or a heavy-flavor jet 
produces leptons via semileptonic decay. Non-VF I^t can 
result from mismeasurements of energy or semileptonic 
decays of heavy-flavor quarks. Since the I^t mismeasure- 
ment is usually not well modeled in detector simulation, 
we directly estimate the contribution of non-VF events 
from the lepton-|-jets data before ^-tagging is applied, 
known as the pretag sample. 

Generally, the bulk of non-W^ events are characterized 
by a nonisolated lepton and small I^t- The lepton isola- 
tion I is defined as the ratio of the sum of the calorimeter 
energy inside a cone of Ai? = 0.4 around the lepton direc- 
tion, excluding the lepton energy, and the lepton energy 
itself. The quantity / is small if the lepton is well-isolated 
from the rest of the event, as typified by a true leptonic W 
decay. This feature is used to extrapolate the expected 
non-W contribution into our signal region, namely, small 
/ and large Ifx- In extracting the non-iy background 
contribution from data, we make the following two as- 
sumptions: lepton isolation and I^t are uncorrelated in 
non-W events, and the 6-tagging rate is not dependent 
on I^T in non-M^ events. The level at which these as- 
sumptions are justified determines the assigned uncer- 
tainty. The contributions from tt and M^-|-jets events are 
subtracted according to the calculated cross sections for 
those processes. 

To validate the method and estimate the relevant sys- 
tematic uncertainties, we vary the boundaries of the sig- 
nal and background regions. The observed deviations 
imply a 25% systematic uncertainty in the non-II^ back- 
ground yield, assigned conservatively for both the pretag 
and tagged estimates 0, ll| . 

The contribution of non-iy background passing the 
NN 6-tagging filter is determined by using data in the 
background region (/ > 0.2 and ^t>20 GeV), which 
has event kinematics similar to non-VF events in the sig- 
nal region because lepton isolation is the only difference 
between the two regions. The non-W^ estimate after ap- 
plying the NN 6-tagging filter is scaled by the ratio of the 



events before NN 6-tagging filter to the events after NN 
6-tagging filter in this background region; this assumes 
the NN 6-tagging filter is uncorrelated with the isolation. 

The non-M^ estimate for events with at least two h- 
tags is obtained by measuring the ratio of the number 
of events with at least one 6-tag to the number with at 
least two 6-tags in the background region and applying 
the ratio to the estimate of tagged non-W^ events in the 
signal region. 

For the plug region, the MET-I-PEM trigger is used, 
which means the above method is not valid, because of 
the I^T trigger bias. The I^t distribution shape differ- 
ence between non-M^ background and other backgrounds 
is therefore used to measure the amount of non-W^ back- 
ground. To model the I^t distribution shape in the non- 
W events, the control samples with electron candidates 
which failed at least two of our standard lepton identi- 
fication criteria are used. We perform a likelihood fit of 
the I^T distribution for observed data using the tiou-W 
and other background shapes [2^ . 



B. W + Mistagged Jets 

The rate at which SECVTX falsely tags light-flavor jets 
is derived from inclusive jet samples by counting the num- 
ber of negative mistags. These are defined as jets tagged 
with an unphysical negative decay length, resulting from 
imperfect tracking resolution. They are assumed to be a 
good estimate of falsely tagged jets, independent to first 
order of heavy-flavor content in the generic jet sample 
and parametrized as a function of 77, number of vertices, 
jet ET^ track multiplicity, z position of primary vertex, 
and total event Et scalar sum. To estimate the posi- 
tive mistag rate, the negative mistag rate is corrected for 
the enhancement of positive tags due to real secondary 
vertices in light-flavor jets and to interactions in the de- 
tector material. This factor is measured in an inclusive 
jet sample by fitting the asymmetry in the vertex mass 
distribution of positive tags over negative tags [1^ . The 
systematic uncertainty on the mistag rate is estimated 
from the differences between observed and expected neg- 
ative tags when applying the same parametrization to 
different jet samples. 

The mistag rate per jet is apphed to W^-|-jets events. 
The total estimate is corrected for the non-VK QCD frac- 
tion and also for the top-quark contributions to the pre- 
tag sample. To estimate the mistag contribution in NN- 
tagged events, we apply the light flavor rejection factor 
of the NN filter 0.35 ±0.05 as measured using light-flavor 
jets from various data and simulated samples. To esti- 
mate the mistag contribution in double tagged events, 
we apply the mistag rate to all untagged jets in ly -I- 1 
6-tagged jet events. This method accounts for the case 
of one real 6-tag plus one mistag and also for the double- 
mistag case in which both 6-tagged jets are not due to 
6-hadron decays. 

For jet probability 6-tagging, the mistag rate is derived 
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from inclusive jet samples as a function of 77. z position 
of the primary vertex, jet E^, track multiplicity, number 
of vertices, and total event scalar Et- The mistag rate 
probabilities are derived from jets with negatively signed 
impact parameters and applied in the same way as for 
the SECVTX algorithm \u\ . 

C. M^+Heavy Flavor 

The Wbb, Wcc, and Wc processes are major back- 
ground sources after the requirement of ^-tagging. Large 
theoretical uncertainties exist for the overall normaliza- 
tion because current Monte Carlo event generators can 
generate VF-|-heavy-flavor events only to leading order. 
Consequently, the rates for these processes are normal- 
ized to data. The contribution from true heavy-flavor 
production in M^-|-jets events is determined from mea- 
surements of the heavy-flavor event fraction in W^-|-jets 
events and the ^-tagging efficiency for those events, as 



explained below. 

The fraction of VF-|-jets events produced with heavy- 
flavor jets has been studied extensively using an ALP- 
GEN -|- PYTHIA Combination of Monte Carlo simula- 
tions [2^, [13] . Calculations of the heavy-flavor fraction 
in ALPGEN have been calibrated using a jet data sam- 
ple, and measurements indicate that a scaling factor of 
1.4±0.4 is necessary to make the heavy-flavor production 
in Monte Carlo match the production in W+1 jet events. 
The final results obtained for heavy-flavor fractions are 
shown in Table [H 

For the tagged M^-|-heavy flavor (HF) background esti- 
mate, the heavy-flavor fractions and tagging rates given 
in Tables HI and [ill are multiplied by the number of pretag 
M^-l-jets candidate events (A^^pietag) in data, after correc- 
tion for the contribution of non-W^ {fnon-w), and other 
background events to the pretag sample. The VF-|-heavy 
flavor background contribution is obtained by the follow- 
ing relation: 



Nw+HF ~ /hF • etag • [A'^prctag ' (1 ^ /non-w) — NtOP — NeWk] , (4) 

I 



where /hf is the heavy-flavor fraction, etag is the tag- 
ging efficiency, TVtop is the expected number of tt and 
single top events, and A^ewk is the expected background 
contribution from WW, WZ, ZZ and Z boson events. 

D. Top and Electroweak Backgrounds 

Production of both single top quarks and top-quark 
pairs contributes to the tagged W^-|-jets sample. Several 
electroweak boson production processes also contribute. 
WW pairs can decay to a lepton, neutrino as missing 
energy, and two jets, one of which may be charm. WZ 
events can decay to the signal ivbb or ii^cc final state. 
Finally, Z — > r+r^ events can have one leptonic r de- 
cay and one hadronic decay. The leptonic r decay gives 
rise to a lepton plus missing transverse energy, while the 
hadronic decay yields a narrow jet of hadrons with a 
nonzero lifetime. 

The normalization of the diboson and top produc- 
tion backgrounds are based on the theoretical cross sec- 
tions listed in Table IIIIl the luminosity, and the accep- 
tance and 6-tagging efficiency derived from Monte Carlo 



events [U, [2^ [2^, [s^ . The acceptance is corrected for 
lepton identification, trigger efficiencies, and the z vertex 
cut. The tagging efficiency is corrected by the 6-tagging 
scale factor. 

VI. HIGGS BOSON SIGNAL ACCEPTANCE 

The PYTHIA Monte Carlo generator is used to gen- 
erate the signal samples and study the kinematic prop- 
erties [13 . Only Higgs boson masses between 110 and 
150GeV/c^ are considered because this is the mass re- 
gion, especially up to 135GeV/c^, for which the decay 
H ^ bb dominates. The number of expected WH — > 
ivbb events N is given by 

= e • y" Cdt- a{pp WH) ■ B{H bb), (5) 

where e, J Cdt, <T{pp WH), and B{H bb) are the 
event detection acceptance, integrated luminosity, pro- 
duction cross section, and branching fraction, respec- 
tively. The production cross section and branching frac- 
tion are calculated to NLO precision [^ . The acceptance 
e is broken down into the following factors: 



Cfctag • (Ctrigger ' Clepton ID ' Ckincmatics ' B{W tv)) , (6) 



where e^o, etrigger, Qepton ID, Cfctag, and ekinomatics are effi- primary vertex, trigger, lepton identification, 6-tagging, 
ciencies defined in sequence to meet the requirements of 
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TABLE I: The heavy-flavor fractions /hf, given in percent, for the W -fjets events where IB, 2B refer to number of taggable 
6-jets in the events, with IC, 2C for charm jets. The results from ALPGEN Monte Carlo have been scaled by the data-derived 
calibration factor of 1.4 ± 0.4. 



Jet Multiplicity 1 jet 2 jets 3 jets > 4 jets 
Wbh (IB) (%) 1.0 ± 0.4 2.0 ± 0.8 3.4 ± 1.4 4.6 ± 2.0 
Wbb (2B) (%) - 1.3 ± 0.6 2.5 ± 1.0 3.1 ± 1.8 
Wcc (IC) (%) 7.7 ± 2.4 12.2 ± 4.5 16.4 ± 5.3 18.6 ± 6.9 
Wcc (2C) (%) - 2.0 ± 0.8 4.6 ± 1.8 8.4 ± 3.4 



TABLE IL The 6-tagging efficiencies etag in percent for the various 6-tagging strategies on individual W-|-heavy-flavor processes. 
Categories IB, 2B refer to number of taggable 6-jets in the events, with similar categories for charm jets. Those numbers include 
the effect of the data-to-Monte Carlo scale factors. 



Jet Multiplicity 


1 jet 2 jets 3 jets 


> 4 jets 




1 SECVTX and NN &-tag (%) 




Wbb (IB) 


27.5 ± 1.5 28.1 ± 1.3 26.7 ± 1.4 


26.9 ± 3.7 


Wbb (2B) 


26.2 ± 1.2 24.1 ± 1.2 


22.6 ± 1.3 


Wcc (IC) 


4.2 ± 0.2 4.6 ± 0.3 4.9 ± 0.3 


5.2 ± 0.7 


Wcc (2C) 


6.3 ± 0.4 6.6 ± 0.4 


6.9 ± 0.6 




> 2 SECVTX b-tag (%) 




Wbb (2B) 


16 ± 2 19 ± 2 


19 ± 3 


Wcc (2C) 


1.0 ± 0.5 2.0 ± 1.0 


2.0 ± 1.0 


1 SECVTX -1- Jet Probability 6-tag (%) 




Wbb (2B) 


10.1 ± 1.2 11.1 ± 1.3 


12.4 ± 1.5 


Wa5 (2C) 


1.6 ± 0.2 2.3 ± 0.3 


3.2 ± 0.4 



TABLE IIL Theoretical cross sections and uncertainties for 
the electroweak and single top backgrounds, along with the 
theoretical cross section for ti at mt = 175GeV/c^ [sj]. The 
cross section of r^r^ is obtained in the dilepton mass 

range m^r > 30GeV/c^ together with a fc-factor (NLO/LO) 
of 1.4 [12 

Background Theoretical Cross Sections 



WW 
WZ 
ZZ 

Single top s-channel 
Single top f-channel 
^ t' 
tt 



12.40 ± 0.25 pb 
3.96 ± 0.06 pb 
1.58 ± 0.05 pb 
0.88 ± 0.11 pb 
1.98 ± 0.25 pb 
265 ± 30.0 pb 



6.7: 



pb 



and event selection criteria. The major sources of inef- 
ficiency are the lepton identification, jet kinematics, and 
6-tagging factors; each has an efficiency between 30 and 
45%. The factor 

fzo is obtained using the vertex distri- 
bution from the minimum bias data, Ctrigger is measured 
using a clean W —f Iv data sample, obtained from differ- 
ent triggers after applying more stringent offline cuts, and 
eiopton ID is calculated using Z U observed data and 



Monte Carlo samples, ettag is measured in a 6-enriched 
sample from semileptonic heavy flavor decay. The total 
signal acceptances are shown in Fig. [T] for the selected 
6-tagging options as a function of Higgs boson mass. In- 
clusion of forward electrons in this analysis increased the 
overall acceptance by 10%. 

The total systematic uncertainty on the acceptance 
stems from the jet energy scale, initial and final state 
radiation (ISR and FSR), lepton identification, trigger 
efficiencies, and 6-tagging scale factor. A 2% uncertainty 
on the lepton identification efficiency is assigned for each 
lepton type (CEM electron, PEM electron, CMUP and 
CMX muon), based on studies of Z boson events. For 
each of the high pT lepton triggers, a 1% uncertainty is 
measured from backup trigger paths or Z boson events. 

The initial and final state radiation systematic un- 
certainties are estimated by changing the parameters in 
the Monte Carlo related to ISR and FSR from nominal 
values to half or double the nominal (sj]. Half of the 
difference from the two samples is taken as the system- 
atic uncertainty. The uncertainty in the incoming par- 
ton energies relies on the eigenvectors provided in the 
PDF fits. A NLO version of the PDFs, CTEQ6M, pro- 
vides a 90% confidence interval of each eigenvector [3a |. 
The nominal PDF value is reweighted to the 90% confi- 
dence level value, and the corresponding reweighted ac- 
ceptance is computed. The differences between nominal 
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FIG. 1: The total acceptance e for the process WH — > £i/bb in 
Vy+2jet bin for the selected 6-tagging strategies as a function 
of Higgs boson mass. 



e, and 

3)[ia 



and reweighted acceptances are added in quadrature, and 
the total is assigned as the systematic uncertainty 

The effect of the jet energy scale uncertainty (JES) 
is calculated by shifting the jet energy scale by ztlcr 
in the WH Monte Carlo samples. The deviation from 
the nominal acceptance is taken as the systematic uncer- 
tainty. The systematic uncertainty on the 6-tagging effi- 
ciency is based on the scale factor uncertainty discussed 
in Sec. HvXl and [TVC] When the NN 6-tagging fiher is 
applied, the scale factor uncertainty is added to that of 
SECVTX in quadrature. The total systematic uncertain- 
ties for the selected 6-tagging options are summarized in 
Table HYl 

The expected number of signal events is estimated by 
Eq. [5]at each Higgs boson mass point. The expectations 
for the selected 6-tagging strategies are shown in Table IVl 



VII. NEURAL NETWORK DISCRIMINANT 

To further improve the signal to background discrim- 
ination after event selection, we employ a neural net- 
work [2^ 1 trained on a variety of kinematic variables to 
distinguish the VF-|-Higgs events from the backgrounds. 

The neural network on the samples of simulated events 
using a mixture of 50% signal and 50% of backgrounds is 
trained. The background composition is chosen to have 
equal amounts of 14^66, tt, and single top, which provides 
a maximum sensitivity over a wide range of input condi- 
tions. 

To optimize the neural network structure, we use an 



iterative procedure to determine the configuration that 
best discriminates signal from the background, and uses 
a minimal number of input discriminants. This is done 
by first determining the best one- variable neural network 
from a list of 76 possible variables, based on the kinematic 
distributions of the two jets, lepton, and I^t in the events 
(including correlations among these objects). The opti- 
mization algorithm keeps this variable as an input and 
then loops over all other variables to determine the best 
two- variable neural network. The best N- variable neural 
network is finally selected once the (NH-l)-variable neu- 
ral network shows less than 0.5 % improvement. The 
criteria for comparing neural networks is the testing er- 
ror defined by how often a neural network with a given 
configuration incorrectly classifies signal and background 
events. This optimization of the neural network struc- 
ture was done for Higgs mass of 120 GeV/c^, and we use 
the same structure to train separate neural networks for 
Higgs masses of 110, 115, 120, 130, 140, and 150 GeV/c^. 
Retraining neural networks with different signal masses 
keeps the neural network improvement almost constant 
as a function of the Higgs mass. 

Our neural network configuration has 6 input variables, 
11 hidden nodes, and 1 output node. The output of the 
neural network has a value from to 1 that will pro- 
vide discrimination between the signal and background 
hypotheses. The signal-like events yield a high neural 
network value and the background like events yield a low 
neural network value. 

The 6 optimal inputs are as follows: 



M. 



33+ ■ 



the invariant mass calculated from the two 



jets in the VF-|-2jets event. Furthermore, if there 
are additional loose jets present, the loose jet that 
is closest to one of the two jets is included in this in- 
variant mass calculation, if the separation between 
that loose jet and one of the jets is Ai? < 0.9. This 
definition gives greater discriminating power for in- 
stances where the selected 6-jet spreads outside the 
cone. 

• ^ E't (Loose Jets): the scalar sum of the transverse 
energies of the loose jets. 

• pt Imbalance: the difference between the scalar 
sum of the transverse momenta of all measured ob- 
jects and the I^t- Specifically, it is calculated as 



Ml^™: the invariant mass of the lepton, I^t, and 
one of the two jets, where the jet is chosen to give 
the minimum invariant mass. For this quantity, the 
Pz component of the neutrino is ignored. 



• AR{lepton-i'rnax)'- the Ai? separation between the 
lepton and the neutrino, where the Pz of the neu- 
trino is taken by choosing the largest \pz\ of solu- 
tions from the quadratic equations for the W mass 
(80.42 GcV/c^) constraint. 
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TABLE IV: Systematic uncertainties on signal acceptance for the selected 6-tagging requirements. 
6-tagging category LeptonID Trigger ISR/FSR JES PDF 6-tagging Total 



ST with NN filter 


~ 2% 


< 


1% 


2.9% 


2.3% 


1.2% 


3.5% 


5.6% 


ST+ST 


~ 2% 


< 


1% 


5.2% 


2.5% 


2.1% 


8.4% 


10.6% 


ST+JP 


~ 2% 


< 


1% 


4.0% 


2.8% 


1.5% 


9.1% 


10.5% 



TABLE V: Expected number of WH Ivhh signal events in central region for the selected b-tagging options. All systematic 
uncertainties are included. 

b-tagging category 110 GeV 115 GeV 120 GeV 130 GeV 140 GeV 150 GeV 

Pretag 9.41±0.61 7.92±0.52 6.35±0.41 3.99±0.26 2.02±0.13 0.78±0.05 

ST with NN filter 2.63±0.22 2.20±0.18 1.70±0.14 1.08±0.09 0.55±0.05 0.21±0.02 

ST+ST 1.18±Q.14 1.00±0.12 0.85±0.10 0.55±0.07 0.27±0.03 O.lOiO.Ol 

ST+JP 0.89+0.11 0.76+0.09 0.60+0.07 0.40+0.05 0.20+0.02 0.08+0.01 



• Pt{W + H): the norm of the vector sum of trans- 
verse momentum of the lepton, I^t and two jets. 

Distributions of all these quantities are checked for 
both the pretag and tagged samples. The simulated back- 
ground events match to the real data well. 



VIII. RESULTS 

A. Counting Results 

The observed number of events in data is compared to 
the expected background in Fig. [21 as a function of jet 
multiplicity. Results are shown for the single and dou- 
ble ^-tagged categories separately. Table |VT] shows the 
composition of + 2 jet data in each fe-tagging cate- 
gory and in the central lepton and plug lepton regions, 
respectively. 

The observed number of events in the data and the 
standard model background expectations are consistent 
in each ^-tagging category. 



B. Limit on Higgs Boson Production Rate 

We apply the neural network discriminant discussed in 
the previous section to the samples of simulated events 
and obtain the distributions of the network output for 
all the processes considered. The expected distributions 
are compared to the data observed in single and double 
^-tagged categories as shown in Fig. [31 We use a binned 
likelihood technique to fit the observed neural network 
output distributions in the three fe-tagging categories to 
test for the presence of a WH signal. 

The number of events in each bin follows the Poisson 
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FIG. 2; Number of events in the central and plug regions 
as a function of jet multiplicity for events with exactly one 
SECVTX &-tag applying the NN 6-tagging filter requirement 
(left) and for events with at least two SECVTX 6-tagged jets 
or one SECVTX b-tagged jet plus one jet probability b-tagged 
jet (right). 
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TABLE VI: Predicted sample composition and observed number of VK + 2 jet events in the central and plug regions with the 
selected 6-tagging options. 



Central region Plug region 



Pretag Events 




32242 






5879 




6-tagging 


ST+ST 


ST+JP 


ST+NN 


ST+ST 


ST+JP 


ST+NN 


Mistag 


3.88±0.35 


11.73±0.92 


107.1±9.38 


1.00±0T8 


3T8+0.49 


28.47+3.30 


Wbb 


37.93±16.92 


31.15±14.03 


215.6±92.34 


7.40±3.96 


6.23+3.37 43.09+12.33 


Wcc 


2.88±1.25 


7.87±3.43 


167.0±62.14 


0.96±0.49 


1.00+0.81 


o3.3/+y.OO 


tt(6.7pb) 


19.05±2.92 


15.56±2.39 


60.68±9.30 


2T4±0.34 


1.79+0.31 


7.17+1.00 


Single top(s-ch) 


6.90±1.00 


5.14±0.75 


14.38±2.09 


0.69±0T0 


0.51+0.08 


1.53+0.20 


Single top(t-ch) 


1.60±0.23 


1.87±0.27 


29.57±4.33 


0.22±0.04 


0.24+0.04 


3.54+0.47 


WW 


0.17±0.02 


0.93±0.11 


15.45±1.91 


O.OliO.Ol 


0.12+0.04 


3.00+0.20 


WZ 


2.41±0.26 


1.84±0.20 


7.59±0.81 


0.58±0.06 


0.42+0.05 


1.62+0.09 


ZZ 


0.06±0.01 


0.08±0.01 


0.31±0.03 


O.OOiO.Ol 


0.01+0.01 


0.02+0.00 


Z ^ TT 


0.25±0.04 


1.29±0.20 


7.27±1.12 


O.OOiO.Ol 


0.01+0.01 


0.24+0.03 


non-W' QCD 


5.50±1.00 


9.55±1.73 


184.7±33.04 


1.16±0.44 


1.51+0.55 


18.34+5.54 


Total Background 


80.6±18.8 


87.0±18.0 


809.6±159.4 


14.2±4.0 


15.5+3.6 


140.4+16.9 


WH signal (120 GeV) 


0.85±0.10 


0.60±0.07 


1.70±0.14 


0.09±0.01 


0.06+0.01 


0.20+0.01 


Observed Events 


83 


90 


805 


11 


13 


138 



distribution 

P^{n,,f,,) = '"' (l = l,2,...,7Vbin), (7) 

where ni,fj,i, and iVbin represent the number of observed 
events in the i-th bin, the expectation in the i-th bin, and 
the total number of bins. The Higgs production hypoth- 
esis is constructed by setting /i^ to fii = Si + bi, where 
Si and bi are the number of signal and expected back- 
ground events in the z-th bin. This quantity Si can also 
be written as a product 

s^ = (t{pP W^H) ■ B{H -^bb)-e- j Cdt- f^" , (8) 

where is the fraction of the total signal which lies 

in the i-th bin. In this case, a{pp W^H) ■ B{H bb) 
is the variable to be extracted from data. The likeli- 
hoods from the three 6-tagging categories are multiplied 
together. The systematic uncertainties associated with 
the pretag acceptance, luminosity, and the 6-tagging ef- 
ficiency scale factor are considered to be fully correlated 
among the three selection categories. Background un- 
certainties on the heavy-flavor fractions are completely 
correlated among VF+HF backgrounds and background 
uncertainties on 6-tagging scale factor are also completely 
correlated among all backgrounds. The systematic uncer- 
tainties associated with the shape of network output arc 
also studied and found to have a negligible impact on the 
final results. 

Since we observe no excess over the background predic- 
tion, an upper limit on the Higgs boson production cross 
section times branching fraction (T{pp W^H) ■ B{H 



bb) is extracted by using a Bayesian procedure. We as- 
sume a uniform prior probability for a ■ B and integrate 
the likelihood over all parameters except for a ■ B. The 
95% confidence level upper limits on tr • S are obtained 
by calculating the 95*'* percentile of the resulting distri- 
butions. 

To evaluate the sensitivity of the analysis, background- 
only pseudocxperiments are used to calculate the ex- 
pected limit in the absence of Higgs boson production. 
Pseudo data are generated by fluctuating the individ- 
ual background estimates within their total uncertainties. 
The expected limit is defined as the median of the 95% 
confidence level upper limits of 1000 pseudocxperiments. 

The observed limits as a function of the Higgs boson 
mass are shown in Fig.[3]and Table IVlTl together with the 
expected limits determined from pseudocxperiments. We 
set 95% confidence level upper limits on the production 
cross section times branching fraction ranging from 1.2 
to 1.1 pb or 7.5 to 101.9 times the standard model expec- 
tation for Higgs boson masses from 110 to 150GeV/c^, 
respectively. 



IX. CONCLUSIONS 

We have presented a search for the standard model 
Higgs boson in the (,vbb final state expected from WH 
production at CDF. The search sensitivity is improved 
significantly with respect to previous searches, by about 
60% more than the expectation from simple luminosity 
scaling. The main improvements are using jet probabil- 
ity 6-tagging, a multivariate neural network technique to 
further enhance sensitivity to the signal, and increasing 
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TABLE VII: Observed and expected upper limits on (j{pp 
WH) ■ B{H bb) at 95 % C.L. compared to the standard 
model production rate calculated at NLO. 

Higgs Mass Upper Limit (pb) Upper Limit/SM 
GeV/c^ Observed Expected Observed Expected 



110 


1.2 


1.2 


7.5 


7.8 


115 


1.2 


1.1 


9.0 


8.7 


120 


1.1 


1.1 


10.2 


10.5 


130 


1.1 


0.9 


17.9 


15.2 


140 


1.2 


0.8 


40.1 


28.7 


150 


1.1 


0.8 


101.9 


70.9 



the acceptance for signal events by including leptons in 
the forward region of the detector. These improvements, 
along with a dataset of 1.9 fb~^, allow us to set a 95% con- 
fidence level upper limit on the production cross section 
times branching fraction that ranges from 1.2 to 1.1 pb 
or 7.5 to 101.9 times the standard model expectation for 
Higgs boson masses spanning from 110 to 150GeV/c^, 
respectively. 

This channel is an important component of the com- 
bined search for the Higgs boson at the Tevatron [36| . 
With the Tevatron expected to deliver a factor of 4 more 
data, it is possible that the combined search with the 



complete dataset will be sensitive to the standard model 
Higgs cross section in this low mass range, which is of 
particular interest according to indications derived from 
standard model fits to electroweak observables H. 
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FIG. 3: Neural Network output distribution in W+2 jets 
events for exactly one SECVTX 6-tagged jet that passes the 
NN 6-tagging filter (left) and events for ST+ST and ST+JP 
double ^-tagging categories (right). The contributions of the 
various background sources from the central plus plug region 
are shown in histograms while the hatched box represents the 
background uncertainty. The expected signal for a 120 GeV 
Higgs boson (multiplied by a factor of 10) is shown by the 
solid line. 
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FIG. 4: The observed (solid) and expected (dashed) 95% con- 
fidence level upper limit on cr(pp WH)-B{H bb) relative 
to the standard model expectations with an integrated lumi- 
nosity of 1.9 fb^^ along with the one and two a bands of the 
distributions of expected outcomes from pseudoexperiments. 



